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are shown in Fig. 6. Note that nearly the 9ame trend (to within 
experimental error) is observed for both powdered and block shale 
samples. Apparently, the degradation process is nearly 
independent of particle size. 

To confirm that this agreement was not just fortuitous, we 
measured the effects of slow heating conditions (2°C/hr) on 
material 3.2-, 7.2-, and 11.8-cmdiam. To within experimental 
error, all particle sizes gave the same yield. A summary of 
yield data from the these experiments is given in Table 2 and 3. 

During previous experiments 5 on powder shale samples, we 
observed similar oil degradation (i.e., yield loss) from prolonged 
isothermal heating under autogenous conditions. However, sweeping 



TABLE 3. a 
Summary of Oil Shale Powder Experiments (22 gal/ton Material 
under Linear Heating Conditions, i.e., dT/dt - Const). 



External 
Heating flow 
Sample rate, ra£e^ Oil yield 

No. 

RP-1 
TC-22 



C/hr cm /min % FA 

"576 0 78.8 

5.0 0 86.1 

TC-18 5.0 0 85.7 

19.3 0 90.8 

63.5 0 96.1 

185.0 0 98. A 



TC-20 
TC-14 
TC-12 

Fischer 
Assay 

TC-28 
TC-52 
TC-34 
TC-56 
TC-38 
TC-48 



720 0 100.0 

5 2.0 86.7 

5 5.2 89.8 

5 7.0 91.9 

5 15.0 92.9 

5 18.6 97.4 

c 19<3 o 91.0 



a 9.914 vtZ organic carbon in raw shale. 

b A zero flow indicates autogenous conditions, 

c Retort vessel inverted. 
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Che sample with an inert gas reduced the yield loss, since the 
liberated oil was removed from the pyrolysis region before being 
degraded . 

The ability of a sweep gas to reduce oil degradation was also 
observed in the non isothermal (linear heating rate) powder 
experiments. Powdered shale samples heated at 5°C/hr showed a 
near linear increase in yield with low gas sweep rate (<18cm 3 / 
min — i.e. < 1 retort void volurae/min.) — see Table 3 and Fig. 7. 
Hence, the effects of gas sweep rate and heating rate closely 
parallel each other. 

These results indicate that a higher heating rate may allow 
the recovery of more oil in both powders and blocks because it 



Material 1 
nst) . 



il yield 
1 FA 

78.3 
86,1 
85.7 
90.8 
96.1 
98.4 

100.0 

86.7 
89.8 
91.9 
92.9 
97.4 
91.0 




< 

i 




Inert gas flow rate - cm /min 

Fig. 7. Oil yield vs external gas flow rate (argon) for powdered 
oil shale samples heated at 5°C/h. A flow rate of 18 cm /min 
corresponds to about one retort void volume per minute at 400°C 
and 1 bar. 
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produces higher self-generated (i.e., autogenous) gas sweep rates 
during retorting. Furthermore, since the oil products are removed 
in the vapor phase, the intra-particle flow rate will be 
approximately proportional to the heating rate. The low 
permeability of oil shale makes it impossible for an external 
sweep gas to significantly penetrate a shale block. Autogenous 
intra-particle gas sweep and oil vaporization processes are the 
only effective means of removing oil from a block. As a result 
the more rapid the heating conditions, the faster oil is removed 
from the block; consequently, less oil degrades. 

We were also interested in the effect of the shale grade on 
yield. To determine this effect, two 17.2-cm diam blocks, one 
14 gal/ton and another 47 gal/ton, were retorted at 2°C/hr and 
compared with experimental data from the approximately 30-gal/ton 
material (Table 2, Fig. 8). To within experimental error, the 
yield is independent of grade over the 14- to 49-gal/ ton range. 

In summary, the experimental data show that oil degradation 
depends principally on heating rate (as it affects the autogenous 
gas sweep rate) and not on particle size or grade for the material 



Grade - gal/ton 
10 20 30 40 50 

90 



i — i — i — [ — i — i — i — i — | — i — i — i — i — | — i — ' — r 

7.0 cm 



17.2 cm 
< |- Mean expt valve 



- 80 

CD 
>- 



\ r 




3.2 cm' 



v, 



8 cm 



t 



70 1 1 1 1 1 I i I — i — i — I — i — i — i— i- 



41.7 83.4 125.1 166.8 208.5 

Grade - 1/t 



Fig. 8. Oil yield (as percent of Fischer Assay) vs grade 

(gal/con, (1/t)), for oil shale particles heated at a rate of 
2 e C/h. 
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used here. (Note that particle size may have an important effect 
on oil degradation because of combustion. Oxidation, however, 
is an extra-particle effect - here we discuss only in trapar t icle 
oil degradation. In a later section, the possible effect of 
heating rate on oil oxidation is discussed.) 

^ntraparticle oil degradation - cracking vs coking reactions 

A realistic picture of oil release from oil shale can be 
developed only when one understands the mechanism of intraparticle 
oil degradation. By determining the mechanisms that lead to major 
oil loss, one can then attempt to design processing conditions 
to minimize loss effects. 

A recent study 5 of the effects of thermal histories on oil 
shale pyrolysis showed that the yield was determined by the time- 
temperature history of the liberated oil and not by the thermal 
history of the kerogen. 5 Using this work as a basis, we have 
confined our investigation to the loss in yield due to reactions 
in the liberated oil. 

The principal nonoxidative means by which the liberated shale 
oil can be degraded are cracking and coking reactions. Because 
coking and cracking often have slightly different meanings, 16 ' 19 
we defined the terms as they will be used here. By cracking 
we mean vapor phase bond fission reactions that eventually lead 
to a distribution of molecular units (mostly units smaller than 
the original molecule) plus, possibly, some carbonaceous residue. 
Coking we define as liquid or condensed phase reactions resulting 
in the fusion of two or more molecular species with the ultimate 
formation of a carbonaceous product plus minor amounts of lower 

molecular weight gases. 

An important difference in the preceding two mechanisms is 
the distribution of final products. As defined, cracking 
reactions lead to greater evolution of lower molecular weight 
hydrocarbon fragments (CH^ t C^, etc.), whereas coking reactions 
produce mainly a carbonaceous residue. 
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Experiments reported in Appendices A and B show that coking 
is the principal intraparticle oil degradation process. Evidence 
for this comes from several observations, principally: 

• Product distribution from powder and block retorting 
experiments (Appendix A). 

• Shale oil distillation at different heating rates 
(Appendix B) . 

Proposed oil production mechanism 

From the foregoing experimental data and that in the 
appendices, it is possible to formulate a simple mechanistic 
picture of oil production and degradation: 

(Kerogen-bitumen) -rrT > ( oil )j + carbonaceous + gas 

/ \ residue m 

V Vl + v 

coke + gas ^ oil produced. 



The term (oil)^ refers to the liberated oil in-place 
(intraparticle). The in-place oil is subjected to two competing 
processes: oil degradation leading to coke and gas formation, 
and oil removal from the particle giving recovered product. In 
the event of negligible degradation the oil produced would equal 
100% of Fischer Assay. 

The experimental results in the preceding section show that 
coke is formed primarly from the oil liquid or condensed phase. 
Therefore, one can rewrite the oil generation and degradation 
steps given in Eq. (I) as follows: 

Oil k d ^ _1 



precursors 



> oil + oil* 

M 

k / \ k. / k (2) 

coke *oil produced. 

+ 
gas 
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Since we are concerned only with reactions involving oil, 
the term oil precursors refers to that portion of the kerogen 
plus bitumen that decomposes to oil (i.e., condensible product 
at 0°C). 

The stoichiometrics (Table A) for the reactions given in 
Eqs. (1) and (2) were determined from experiment. Details 
leading to this formulation are given in Appendix C. 

In the proposed mechanism, the oil precursors decompose at 
some rate k giving a distribution of oil in the gas (oil 8 ) and 
liquid (oil 1 ) phases. This distribution of products is some 
function of temperature (X(T)) times the total oil present. 



TABLE 4. 

Stoichiometrics for the Reactions given in Eqs. (1) and (2) 

of the Text. 



Oil generation 

CK 1.54 N 0.028 S O.0O5°O.O46 
Raw shale 



Oil degradation 
CH 1. 56 N 0. 020*0. 002 



-> 0.2 CH n , c S n n . 0 Carbonaceous 

U.UOO U.UlcS 

residue 



+ °' 725 CH 1.56 N 0.020 S 0.O02 011 



+ 0.018 C0 2 
0.019 H 2 
0.013 CH 4 
0.014 C 3 _ 2 H 7-5 



+ 0.010 H 2 0 



> Gas 



Water 



* °' 83 ^0. 2*0. 024*0.0024 C ° ke 



+0.17 CH^ 
+ 0.36 m 



Gas 



Oil removal 
CH 1.56 N 0.020 S 0.002 



k, 1. 56*0. 020 5 0. 002 
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oil g - (X(T)) oil (3a) 
oil 1 - (1 - X(T)) oil (3b) 

The term oil 8 here refers to the material that is above its 
boiling point, i.e., material that cannot coexist in liquid-vapor 
equilibrium at 1 bar and the retort temperature, T. The oil 
vapor in equilibrium with a liquid oil phase is considered 
separately in the k^ rate step. 

Since the void space in the sample is small, the vapor 
products are assumed to exit the material at some rapid rate 
(k ) giving "produced" oil. 

Inside the shale material, the oil in the liquid phase is 

subject to two competing reactions: rate, k , which leads to 

c 

coke formation, and k^ which leads to oil produced. The rate of 
coking, k , is for a reaction and hence varies as exp(-E/RT), 
whereas k^ is essentially independent of temperature, depending 
only on the gas flow rate (or heating rate). 

The process governed by the rate k^ can be considered an 
evaporation step and represented by: 

oil 1 t oil vap + oil produced, (4) 
1 2 

where step 2 depends on the flow rate. Oil vap is the vapor 
material in coexistence with some liquid phase, and hence is 
distinguished from oil®, the oil totally in the gas phase — i.e., 
above its boiling point. 

From the above proposed mechanism it is apparent that 
increasing the flow rate (by either rapid heating or, in the case 
of powders, introducing an inert gas) increases the rate of oil 
production and reduces the likelihood of competing coking 
reactions . 

To carry out meaningful calculations using the mechanism in 
Eq. (2) , one must know the values of the various rate constants 
as well as the boiling point distribution for the produced oil 
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[i.e., X(T)]. Appendices D and E give detailed descriptions of 
how values were determined for these various quantities. The 
rate constants are summarized in Table 5, 

Although the mechanism is over-simplified, the complexity of 
the system precludes a more detailed picture with the data now 
available. Nevertheless, this simple mechanism helps immensely 
in forming a picture of the processes that may be occurring, as 
well as in defining retorting conditions that may lead to low 
yields . 

Analysis of block retorting by a simple mathematical model 

In this section we analyze the results from the block 
experiments using a simple mathematical model. The model 
incorporates the kinetic results from the powder experiments 
[Eq. (2) and Tables 4 and 5] and thermal data reported in the 
literature. The model calculates the volume of evolved oil, the 
degree of intraparticle oil degration, the coke distribution, 
and the thermal profile through the material. By carrying out 
such a numerical analysis we can check the ability of our 
simplified model to describe the experimental results. If there 
is good agreement, the model can be used to gain further valuable 
insight on processing conditions and block sizes leading to 
optimal oil recovery. 



TABLE 5, 

Summary of Rate Constancs for Proposed Mechanism 

[Eq. (2) of text]. 
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k A (s 1 ) = 2.8 x io 13 ' 

" -1 7 

k (s ) = 3.1 * 10 ■ 

k^(s ) = a # b 



exp(-52400/RT) 
exp(-350OO/RT) 



'Cf 1 



k >>> k. , k 
g 1' c 



a = 0.12 
b = heating rate 



Activation energies are in cal/raole. 
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Particle size 

To make the problem numerically tractable, we assume the 

block car be treated as a sphere (i.e., one dimensional). It 
20 

has been shown that particles with equal dimensions in three 
directions can be analyzed by assuming that the particle is a 
sphere with a radius given by 

3V 

r s " AT 2 • (5) 

P 

where and A p are the values for the volume and area of the 
particle, respectively. The 17.2-cm diam blocks have a length- 
to-diameter ratio of about one and, therefore, can be treated 
using this assumption. 



Governing equations 

The equation describing energy conservation for a solid 
spherical particle with a chemical heat source is given in 
Part A of Table 6. The three terms on the right refer to energy 
transport by conduction, convection, and chemical reaction, 
respectively . 

In the analysis given here the convective heat transport 
term is assumed to be small and hence is neglected. In 
Appendix F this assumption is shown to be quite valid for the 
heating rates and grade of material used here. 

The organic reactions and rate of heat generation during 
retorting of an oil shale particle (below the dolomite and 
calcite decomposition temperatures) can be described by the set 
of rate equations given in Parts B and C of Table 6, These 
equations come directly from the experimentally derived mechanism 
for oil generation and degradation given in the previous section 
[Eq. (2)] and the stoichiome tries in Table A. 

The density of the shale block at a particular temperature 
is given by the expression 
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p s (T) ' °min + °org ' O 

Assuming negligible decomposition of mineral species (T < 600°C), 
the block density can be written as: 

P s (T) = °°s ~ °vol (T) <*> 

where p° is the initial block density. The volatile mass loss 
under zero degradation conditions (i.e., 100% FA, by definition) 
is given by P yQ ^ where 

P vol (T) « 0.78 (1 - F) p° b . ( 8 ) 

The coefficient 0.78 represents the mass fraction of volatile 

2 1 

matter released per unit weight of initial organic material. 

At the end of retorting the intraparticle oil has either vaporized 

or coked and hence: 

f o 

P = P - P , (9) 

s s vol 

This system of equations describing the particle retorting 

process was solved using a partial differential equation routine 

22 

implemented by Sincovec and Madsen. The particle was zoned 
into 21 equal-volume segments. The assumed boundary conditions 
are 



|f) =0 (10) 



and 



T l = T Q + bt. (11) 



s 

The second boundary condition, Eq . (11), is set by the experiment, 
where T Q is the initial block temperature and b is the heating 
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experiment , 
heating 




rate determined from the thermocouple readings near the block 

edge (r/r » 0.97) . 
s 

The other necessary initial conditions include the block 
density and concentration of organic matter. The block density 
was determined experimentally (see experimental section), and 
the concentration of organic matter was calculated from the oil 
grade using the correlation reported by Stanfield et al. 2 ^ 

Oil shale physical properties 

The heat capacity of the shale is taken to be: 



C = 



[(216.8 + 120.9 W )-(1.0 - F) 



+ 197.8 F] + [(0.1478 
+ 1.329 W )•(!., 0 -F) 
+ (0.220 F)]-[T - 298] 



(12) 



where is the weight fraction of kerogen in the raw shale and F 
is the fractional conversion (decomposition) of kerogen at 

temperature T. This equation is also from the analysis by 

24 25 
Carley of the Wise et al . data. 

The thermal conductivity of the shale is assumed to be given 



bv 



k = 1.25 * 10" A -(1 - F) 



-4 



+ 0.83 * 10 F , 
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(13) 



based on data reported by Tihen et al. for 30-gal/ton material. 

Model calculation vs experimental results 

It is important to note that there are no free parameters in 
che model. All input to the model was from experimental results 
given here (e.g., kinetics from powder shale experiments) or data 
previously reported in the literature (e.g., heat capacity, 
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thermal conductivity, etc.). Therefore, a comparison of the 
calculated and experimental results provides a sensitive test of 
the model. 

Oil evolution and temperature distribution 

Plotted in Figs. 9 and 10 are calculated and experimental 
results showing the oil evolution, total oil yield, and 
temperature data for a set of block retorting experiments. The 




Block surface temperature - C 

Fig. 9. Oil evolution vs block surface temperature calculated 
using mathematical model and observed experimentally for surface 
heating rates of (a) 2, (b) 6, (c) 18, and (d) 60°C/h. Blocks 
are 17.2-cm diam and 29 to 33 gal/ton. The experimental points 
are from Figure 5. 
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Fig. 10, Temperature difference between block surface and center 
vs block surface temperature, calculated using mathematical 
model and observed experimentally (•) for heating rates of (a) 
2, (b) 6, (c) 18, and (d) 60°C/h. Blocks are 17.2-cm diam. 
The blocks are the same as those referred to in Fig. 9. 
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experimental data are for the 17.2-cm diam blocks heated at rates 
of 2, 6, 18, and 60°C/min [corresponding to sample Nos. 111-38, 
18, 30, and 63, respectively (Table 2)]. The grade of the blocks 
ranged from about 29 to 33 gal/ ton. 

The calculated evolution of oil vs block surface temperature 
agrees well with the experimental results (Fig. 9). The biggest 
discrepancy between the model calculation and experimental data 
is for the high heating rate experiment (l°C/min). The model 
predicts oil release at a slightly higher temperature (-20°C) 
than is observed experimentally. In all other cases, however, 
the experimental and calculated evolution profiles agree to 
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within 10 4 C. Furthermore, the calculated effect of heating rate 
on total oil yield agrees with the experimental results. 

The measured and calculated temperature difference between 
the block surface and center is shown in Fig. 10. The data are 
Plotted vs the block solid surface temperature. The conditions 
correspond to those in Fig. 9. Temperature difference was chosen 
to represent the thermal data, since it provides a more sensitive 
comparison between calculated and experimental results than do 
temperature profiles. 

The agreement between the calculated and observed temperature 
values is quite good. In all cases the model predicts the block 
center temperature to within just a few degrees. 

The biggest discrepancies appear near the end of oil 
evolution; this may indicate the occurrence of other endotherraic 
reactions (e.g., pyrolytic gas evolution from the carbonaceous 
residue or carbonate decomposition) and/or an incorrect value for 
the thermal conductivity of the spent shale. 

Accurate data for the difference between block center and 
surface temperatures were not available for the heating rate of 
2°C/hr. The error in the thermocouple readings was too large to 
quantitatively determine a difference except to note that it was 
less than 10 to 15°C, hence only calculational results are shown 
in Fig. 10(a). 



Carbon distribution in spent shale 

The radial distribution of organic carbon in the spent shale 
blocks was determined at the completion of several retorting 
experiments. Analyses were performed on drill corings taken at 
different radial positions; the procedure used in the coring 
operations is similar to that reported by Mallon and Braun. 27 

Shown in Fig. 11(a) and (b) are typical results for two 
17.2-cm diam cylindrical blocks heated at rates of 2 and 
60°C/hr. The circles represent experimental data, whereas the 
lines are results from model calculations. The interesting 
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Fig. 11. Weight percent organic carbon remaining in spent shale 
vs radial position in a 17.2 cm diameter shale particle. The 
straight line is from a model calculation and open circles 
represent experimental data for 121.8 and 156. A 1/t (29.2 and 
37.5 gal/ton) blocks heated at (a) 2 and (b) 60°C/h, respectively 



feature is that the concentration of organic carbon residue is 
independent of radial position in the block. This phenomena was 
also observed in blocks having a smaller diameter. See, for 
example, Fig. 12 for a 7.2-cm diam block. (It was not possible 
to use the 1-D model for the 7.2-cm block since the length- 
diameter ratio was 2.5. Therefore, the solid line in Fig. 12 
was determined from the observed yield loss and the calculated 
organic carbon converted to char - see Table A2 of Appendix A). 
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Fig. 12. Weight percent organic carbon remaining in spent shale 
vs radial position in a 156.0 1/t (37.4 gal/ton), 7.2 cm diameter 
cylindrical block. The solid line was calculated and the circles 
represent experimentally observed values. The heating rate was 
2°C/h. 



The spatial independence of organic carbon deposition 
indicates that the oil that degrades is apparently not diffusing 
through the block, but instead remains near the site of initial 
formation. This observation supports the proposed mechanism for 
oil degradation [eq. (2)]. The mechanism assumes that the oil 
degradation due to coking occurs in the relatively immobile 
condensed or liquid phase and not in the mobile gas phase; hence, 
the mechanism predicts a spatially independent deposition of 
carbon during retorting (for a near constant heating rate 
throughout the particle). 

From the preceding data it becomes clearer why particle size 
has little effect on oil yield. Since the degradation process 
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is not occurring in the gas phase, longer int raparticle residence 
times (i.e., larger blocks) have little effect on the degradation 
process. This conclusion pertains only to block sizes and 
heating rates similar to those used in this study. One can 
easily think of conditions (e.g., very fast heating rates and/or 
large blocks) where the surface temperature of the shale becomes 
very hot (>700'C) while the inside is still retorting. Under 
these conditions, cracking of the gas phase products diffusing 
from the interior of the block may be an important oil degradation 
process. It is also probable that extrapar t icle oil combustion 
will become a major oil loss mechanism. 

Water evolution during block retorting 

The amount of water released during retorting of the 17.2-diam 
blocks is given in Table 7. The weight of water released is 
about 1.9Z of the total weight loss (to -500°C) from the 29- to 
37-gal/ton material. 

The dynamics of water evolution above 240°C were recorded 
during two of the slow retorting experiments. The water is 
released in approximately equal amounts in two steps (Fig. 13) 



tion 

t diffusing 
of initial 
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it the oil 

mobile 
hase ; hence, 
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rate 

article size 
m process 




Run No, 



III-2 

111-18 

111-30 

111-38 

111-63 

111-55 

111-81 



TABLE 7. 

Quantity of Water Released during Retorting of 
17.2-cm diam Blocks 



Grade, 
gal/ton 



H 2 0 yield 
ml H 9 6/kg raw shale 



29.0 
32 .4 
33.3 
29.2 
37.5 
13.8 
48.8 



25 
21 
17 
17 
17 
12 
14 £ 



Accurate determination of the total water released was 
complicated by the large quantity of oil evolved. Hence, this 
number represents a lower limit on the total water evolved. 
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Fig, 13. Volume of water released (as % of total) vs block 
surface temperature for a 133.5 and 121.0 1/t (32 and 29 gal/ 
ton) shale blocks (17.2 cm diameter) heated at 2 and 6°C/h, 
respectively . 

one below -~240 °C and the second during retorting (i.e., above 
240°C). 

Note the similarities in the evolution profiles for the water 
with those for oil (Fig. 5). The water released at high 
temperatures (>240°C) apparently forms during the kerogen 
decomposition process. The water released below 240°C may come 
from hydrated mineral species or nahcolite in the shale. 

Practical implications of the results 

The results of this study have a great deal of practical 
significance for proposed methods of in- situ shale processing. 
As mentioned in the introduction, there are two principal means 
of processing shale: combustion and hot inert gas retorting. 
The movement of the pyrolysis front during combustion retorting 
is quite fast, the shale heating rates being of the order of 
several degrees per minute. On the other hand, heating rates 
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during hot inert gas retorting may be much 3 lower - only a few 
degrees (or few tenths of a degree) per hour. 

From the results of this study it is apparent that significant 
intraparticle degradation may occur during hot gas retorting, 
whereas relatively little will probably result from combustion 
processing . 

A least-squares fit of the data for oil yield vs logarithm 
of the heating rate (data are from both powders and blocks - 
Fig. 6 and Tables 2 and 3) gives a valuable empirical relation 
that can be used to estimate the intraparticle yield loss at 
various heating rates: 
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YL = 100 - [81.1 + 8.24 (log b) ] , 

b < 190°C/hr 
YL ; 0.0, b > 190°C/hr. 



(14) 



Furthermore, the percent organic carbon converted to char within 
the block can be estimated as 



0C char = -°' 606 



- 19.9, 



(81.1 + 8.24 log b) 
+ 80.54, 
b < 190°C/hr 
b > 190'C/hr. 



(15) 



The above expressions are valid over the range of heating rates, 
particle size, and grade of material used in this study. Caution 
should be used in applying these relationships to very large shale 
particles, since the material used in this study was less than 
18-cm diam. LLL is currently conducting experiments to check' the 
validity of this work when applied to particles approaching 60-cm 
(2 ft) diam. 

The data here also only apply to intraparticle degradation. 
Oxidation of the oil once it leaves the blocks (extraparticle) 
can become important for large size material processed at a high 
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heating rate. Under this condition, a large thermal gradient 
develops, resulting in kerogen pyrolysis occurring on the interior 
of the block while the exterior is being oxidized (i.e., burned). 
LLL is currently developing a detailed retorting model that cakes 
into account both of these processes. 

SUMMARY AND CONCLUSIONS 

The principal conclusion of this study is that particle size 
(up to 17.2 cm) and grade (14 to 49 gal/ ton) have little effect 
on the degree of intraparticle oil degradation during shale 
retorting. On the other hand, the rate at which the shale is 
heated has a pronounced effect on yield. For example, retorting 
at a rate of 180°C/hr gives a yield of -99%, while retorting at 
2°C/hr gives only about 83% yield. 

The principal degradation product is a carbonaceous residue, 
which we term coke. Evidence shows that oil degradation occurs 
mainly in the liquid or condensed phase (coking) and not in the 
vapor phase (cracking). Furthermore, the degree of degradation 
is controlled by the rate of autogenous (i.e., self-generated) 
internal gas sweep through the particle, which in turn is 
controlled by the rate of particle heating. Therefore, higher 
heating rates generate greater intraparticle gas sweep rates 
and less oil degradation. 

A mechanism is proposed to explain the observed degradation 
of oil during retorting. Rate equations characterizing the 
various mechanism steps are determined from experimental data on 
powdered shale and incorporated into a mathematical model 
describing oil shale block retorting. The model contains no free 
parameters. Calculations carried out with the use of the model 
agree very well with experimental results. 
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FOOTNOTES 

'Reference to a company or product name does not imply 
approval or recommendation of the product by the University of 
California or the U.S. Department of Energy (DOE) to the 
exclusion of others that may be suitable. 

SYMBOL LIST 



2 

= particle surface area (m ) 
= shale heat capacity (kcal/kg K) 

= heat capacity of evolved gaseous products (kcal/kg K) 
= activation energy (kcal/mole) 
= dimensionless fraction of decomposed kerogen 
= enthalpy of kerogen decomposition (kcal/kgm) 
char = P ercent organic carbon converted to char 



2 



C 

g 

E 
F 

AH 
OC 



R 



R 



ol - rate of liquid oil formation (intrapart icle) (kg/s-m 3 ) 



>f 

'tain 



R 



R 

op 

T 
T 

o 

V 

P 

w, 

K 

X 

1-X 



VT 



rate of gas phase oil formation (intraparticle) 
(kg/s-m 3 ) 

rate of intraparticle coke formation (kg/s-m 3 ) 
rate of non-condensible gas formation during coking 
(kg/s-m 3 ) 

rate of oil release from the particle (kg/s-m 3 ) 
temperature (K) 

initial particle temperature (K) 

3 

particle volume (m ) 

weight fraction kerogen in raw shale 

fraction of liberated oil totally in the gas phase 

fraction of liberated oil as liquid or in liquid-vap 

equilibrium. 

precent oil yield loss 
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b - heating race (K/s) 

f Q - fractional conversion of kerogen-b itumen to oil 

f =» fractional conversion of degraded oil to coke 

c 

k =■ shale thermoconductivi ty (kcal/s »m* K) 

k, - rate constant for oil generation (s ^) 

d 

k. - rate constant for oil production via liquid 

-1 

vaporization (s ) 
k = rate constant for oil production as initial gas products 

(s l ) 

k = rate constant for coke formation (s ^) 

• c 3 

q = rate of heat generation (kcal/s # m ) 

r = radial position (m) 

r = block radius (m) 
o 

r = radius of equivalent spherical particle (m) 
s 

t - time (s) 

3 

p = density of shale particle (kg/m ) 



3 

P . = mass concentration of mineral matter (kg/m ) 



s 

min 

3 

p = mass concentration of organic matter (kg/m ) 

org 

p° = initial shale particle density (kg/m^) 

s 

f 3 
p = final shale particle density (kg/m ) 

s 

3 

0 = mass concentration of intraparticle liquid oil (kg/m ) 

lo 

p° = initial mass concentration of kerogen + bitumen (kg/m ) 

kb 

p , » mass concentraion of kerogen + bitumen converted to 
vol ^ 

volatile products (kg/m ) 



p v = average mass flux rate of gas products from shale 
g g 2 

(kg/m sec) 
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APPENDIX A. EVIDENCE FOR COKING: ORGANIC CARBON DIST RIBUTION 
IN RETORTING PRODUCTS ~ " 

The amount of organic carbon ' converted to char under various 

retorting conditions is given in Table Al for powders and Table A2 

for four of the block experiments. Stout et al. 5 reported that 

83% of the carbon lost from oil by degradation processes appears 

as char and 17% as gas (mainly CH^). Furthermore, they showed 

that the percent organic carbon converted to char (0C ) could 

char 

be correlated with oil yield by the empirical relation: 



0C 



char - -0.606 • 0Y + 80.514, 



where 0Y is the oil yield as percent FA. Using this correlation, 
a predicted value of organic carbon converted to char has been 
calculated for the experiments in Table Al and A2 . As shown in 
the tables, the agreement between the predicted and observed 
values is quite good. 
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TABLE A2. 

Organic Carbon Converted to Char for Several Block 
Retorting Experiments. 



Sample 


Grade, 
gal/ton 


Heating 
rate 
°C/hr 


Yield, 
X FA 


% organic 
converted 

Measured 


carbon 
to char 

Predicted 5 


III-2 


29.0 


180 


99.0 


20.3 


20.5 


111-38 


29.2 


2 


83.8 


29.3 


29.7 


111-63 


37.5 


60 


96.0 


20.2 


22.3 


111-85 


37.4 


2 


84.9 


29.4 


29.1 


a Block 


dimensions 


are given in 


Table 2 


of main text. 





b Based on work by Stout et al. 5 , organic carbon converted to 
char = (-0.606 x oil yield) + 80.54. 
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The important point is that over 80% of the carbon lost from 
the oil goes to form carbonaceous residue — this strongly indicates 
coking as the major oil degradation process. Furthermore, Stout 
et al. found that the H 2 concentration in the evolved gas 
increased by roughly twofold with a 20% loss in FA yield. There 
was no observable increase in C , n > 1, gas phase hydrocarbons. 
Hence, it is probable that for the experiments reported here, 
little loss in yield results from cracking of the produced oil. 



APPENDIX B. EVIDENCE FOR COKING: SHALE OIL DISTILLATION 
EXPERIMENTS " ~ ' 

Coking of shale oil is readily observed when distilling it 
at various heating rates. Shown in Fig. Bl is the weight precent 
coke produced vs log of the heating rate during distillation of 

shale oil at 1 bar. Also plotted ( ) is the equivalent amount 

of carbonaceous residue formed by oil degradation during retorting 
of shale powders. The dashed curve was obtained using the 
expressions : 

weight % coke = YL • (0.83) • (13.0/13.5) , 
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Heating rate - C/h 



Fig. Bl. (-o-) weight percent coke produced during distillation 

of ^ shale oil at different heating rates and ( ) calculated 

wt% coke formed by retorting powdered shale. 



where YL is the percent yield loss observed at a given heating 
rate (Fig. 6), 0.83 is the fractional conversion of organic carbon 
in oil to organic carbon in coke, and 13.0 and 13.5 are approximate 
empirical weights for 500°C carbonaceous residue (CH) and oil 
(CH^ respectively. 

The agreement between the weight of carbon residue formed 

by retorting powdered shale and coking shale oil is quite good. 

This tends to further support coking as the major Intra-particle 

oil degradation process. 

APPENDIX C. REACTION STOICHIOMETRIES 

We have recently reported the organic carbon distribution 
and elemental composition for the oil, char, and gas produced 
during pyrolysis of oil shale. 5 For a typical Fischer Assay 
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heating schedule (12'C/min), 202 of the organic carbon orginally 
present in the kerogen is converted to char, 72. 52 to oil, and 
7.52 to gas. The oil and char elemental compositions are 

CH 1.56 N O.02O S O.0O3 and ^0 .066 S 0. 022 9 re3 P*ctively and the gas 

is 24.21 CO 26.12 H 18.42 CH. , and 25.92 Z. C by volume 
* *• n>^ n y 

(average composition for E >2 C n - C 3 2 H ? ). On the basis of 
these results and the elemental composition of kerogen,^ 8 one 
can write the following stoichiometric equation for the 
decomposition of kerogen (Fischer Assay yield) : 
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CH 1.5A N 0. 028 S . 005°. 046 



0.725 CH, e£ N rt „S 

1.56 .020 .002 

+ 0.2 CHN A ,,S 

.066 .018 

+ 0.018 CCL 



+ 0.019 H 2 
+ 0.013 CH, 



+0.014 C 32 H 75 
+ 0.010 H 2 0 



Oil 
Residue 



Gas 



Water, 



Only very slight adjustments needed to be made in the observed 

elemental and gas compositions for the equation to balance. The 

water product was added to close the hydrogen balance. 

We also found 5 that during oil degradation 83% of the carbon 

lost from oil appeared as coke and 172 as gas in the form of CH . 

4 

Evolution of considerable quantities of hydrogen was also 
observed. 

The coke formed during oil degradation is lower in hydrogen 
content than the carbonaceous residue from a Fischer Assay. The 
carbonaceous residue has an H/C ratio of about one, whereas the coke 
plus carbonaceous residue from a 792 oil yield experiment has an 
H/C ratio of only 0.70. From the results in Table Al of Appendix A 
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one can calculate a coke H/C ratio of about 0.2. Hence, the oil 
degradation reaction can be written as: 

k 

CH 1.56 N 0.020 S 0.002 ~ > 0-83 CH^ . 2 N Q. 020^0 .002 

+ 0.17 CH^ + 0.36 H 2 

APPENDIX D. RATE CONSTANTS FOR OIL GENERATION AND 
DEGRADATION CALCULATIONS 

The rate constants for the proposed mechanism were determined 
using data from powder experiments. All reactions were treated 
as first order. 

The rate of oil generation has recently been reported by 

29 

Campbell et at. and found to have a rate constant: 

-1 13 
k d (s ) = 2.81 x io exp(-52400/RT). (ID) 

As a first approximation the rate constant for coking, k 

c 

was assumed to be equal to that reported by Raley and Braun fox 
30 

asphalts. The temperature dependence (i.e., activation energy) 
of this rate data, however, is too high and predicts <10% 
difference in yield over a change in heating rate of 720 to 
l°C/hr. The observed yield loss is ~~20% (Fig. Dl). By using a 
lower temperature dependence for the coking rate (E = 35 kcal/mole) 
one can predict the observed yield losses with decrease in heating 
rate. The value of the rate constant at 350°C, however, is still 
assumed to be given by that for asphalt coking (i.e., 1.6 * 10 

sec k is then: 

c 



-5 



k c (a" 1 ) = 3.1 x 10 7 . exp(^|00). (2D ) 

Since k g is much greater than the other rate constants, it 

has no influence on the calculations and hence is arbitrarily set 

equal to a large value such that k » k , , k . k„ . 

g d' c' 1 
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Fig. Dl. Calculated ( ) and observed (•) oil yield from 

powdered 22 gal/ton oil shale vs heating rate. 



Unfortunately, no quantitative information is available for 
k r only the observation that ^ a flow rate a heating rate; the 
proportionality constants are unknown. For a rate of coke 
production given by Eq. (2D) it was found that the value of k 



k » a • heating rate, 
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where 



a - 0.12 °K 



-1 



(3D) 



gives excellent agreement with the observed yield for retorting 
of powders at a heating rate of 60°C/hr. 

Using the above rate constants, the oil boiling point data 
in Appendix E and the mechanism shown in the main text [Eq. (2)] t 
one can calculate the dependence of oil yield on heating rate for 
powder retorting experiments (Fig. Dl). The calculations were 

carried out using a standard ordinary differential equation 

31 

computer routine. 



APPENDIX E. BOILING POINT RANGE OF OIL PRODU CTS BY SIMULATED 
DISTILLATION ~ " " 

The distribution of products in the liquid and vapor phase 
was estimated via a gas chromatographic simulated distillation 
of a Fischer Assay oil (Fig. El). In the numerical calculations 
it is assumed that the weight percent of the material completely 
in the gas phase (i.e., above its boiling point) is given by the 
expression (Fig. El) 

325) • 0.0035, (IE) 
for T > 217°C 

for T > 500°C . 

The material not totally in the vapor phase is designated 
as oil 1 and is assumed to be at some equilibrium vapor-liquid 
distribution (see discussion in main text). 

The rather poor agreement between Eq. (IE) and the distillation 
curve below 300°C (Fig. El) has no effect on the calculations 
since retorting begins above 300°C. 



X(T) ■ 0.38 + (T - 
= 1.00, 
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Fig. El. Simulated distillation curve for Fischer Assay oil. 
Dashed line represents linear approximation of boiling point 
curve used in numerical calculations. 
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APPENDIX F. I NT RAP ARTICLE CONVECTIVE HEAT TRANSPORT 

An accurate account of convective heat transport requires a 
detailed model description of fluid flow through the porous shale 
media. This not only greatly increases the complexity of the model 
but also necessitates the use of several free parameters, since 
little quantitative information is available on the properties 
(porosity and permeability) of the shale. In this appendix it 
is shown that for the case of a constant mass flux rate (similar 
to that encountered during retorting) the effect of convective 
heat transport is small compared to conduction and hence can be 
neglected. 

From the data in Fig. 5 and Table 2 (main text) it can be 
shown that average mass flux rate (p~v~) during retorting of a 

5 a 

17.2-cm diam 37 gal/ton shale block at a heating rate of 60°C/hr 
is approximately 8.7 x i<f 4 kg/m 2 s. Assuming the average heat 
capacity of evolved product is 0.80 kcal/kg °C (that for n-decane) , 
one can then determine the temperature distribution in the block 
with both convective and conductive heat transport. The method 
of solution is the same as that discussed in the main text. In 
the calculations the mass flux at a particular radial position 
in the block is assumed to begin with the onset of kerogen 
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Fig. Fl. Calculated temperature difference between block surface 
and center vs surface temperature. The dashed line is for heat 
transport by conduction and convection, whereas the solid line 
is for conduction only. The surface heating rate is 60°C/hr 
and the particle is 17.2 cm in diameter and 37 g/ton. 
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rig. F2 . Oil evolution vs block surface temperature calculated 
for the cases of heat transport by conduction only (-) and by 

conduction plus convective ( ). The surface heating rate is 

60°C/hr and the particle is 17.2-cm diam and 37 gal/ton. 



decomposition and continue during the entire heating process. 

The results of the above calculations are given in Fig. Fl 
and F2 . Plotted is the temperature difference between the block 
surface and center (Fig. Fl) and amount of oil evolved (Fig. F2) 
vs block surface temperature. Also shown are the results from 
calculations including only conduction. The results show that 
the convective effects are minor compared to conduction and 
consequently can be neglected. 

NOTICE 
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